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^ ' Abstract 

CN ■ 

In this article, we take the point of view that the pentaquark state © + (1540) 
has negative parity, and choose the diquark-triquark type interpolating cur- 
rent to calculate the strong coupling constant qqnk in the QCD sum rules 
£T) I approach. Our numerical results indicate the values of the strong coupling 

constant g<sNK ar e very small, Iffev/d = 0.175 ± 0.084, and the width Tq < 
4MeV, which can explain the narrow width T < WMeV naturally. 
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1 Introduction 

^. 

Several collaborations have reported the observation of the new baryon + (154O) 
Ph! with positive strangeness and minimal quark contents ududs [I]. The existence 

predicted in the chiral quark soliton model, where the + (154O) is a member of 
the baryon antidecuplet 10 The discovery has opened a new field of strong 
interactions and provides a new opportunity for a deeper understanding of the low 
energy QCD. Intense theoretical investigations have been motivated to clarify the 
quantum numbers and to understand the under-structures of the pentaquark state 
+ (154O) The zero of the third component of isospin 7 3 = and the absence of 
isospin partners suggest that the baryon + (154O) is an isosinglet, while the spin and 
parity have not been experimentally determined yet and no consensus has ever been 
reached on the theoretical side. The extremely narrow width below lOMeU puts 
forward a serious challenge to all theoretical models, in the conventional uncorrelated 
quark models the expected width is of the order of several hundred MeV, since the 
strong decay + — > K + N is Okubo-Zweig-Iizuka (OZI) super- allowed. 

In this article, we take the point of view that the quantum numbers of the 
pentaquark state + (154O) are J p = | ,7 = 0,5' = +1, and study its decay 
width within the framework of the QCD sum rules approach jHEHH]. 

The article is arranged as follows: we derive the QCD sum rules for the strong 
coupling constant of the pentaquark state + (154O) Qqnk in section II; in section 
III, numerical results; section IV is reserved for conclusion. 
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2 QCD sum rules for the coupling constant g®NK 



In the following, we write down the three-point correlation function U\ , 

r(p,q) = Jd'xd'ye^e-^iOlTiM^JKiy^emiO), (1) 

where the tjn, jx and rj B = 77q7 are the interpolating currents for the neutron, K 
meson and pentaquark state + (154O) respectively, 

Ve(0) = -±=e abc {<(0)C T5 4(0)} K(0)s e (0)i 75 4(0) - d e (0)s e (0)i l5 u c (0)} ,(2) 

3k{v) = s(y)i~/ 5 u(y) , (3) 
rj N (x) = e^idl^C^d^j^u^x). (4) 

For the pentaquark state + (154O), we use the diquark-triquark type interpolating 
current which can give satisfactory mass and stable magnetic moment [7J |H] . The 
pseudoscalar mesons ir and K can be taken as both Goldstone bosons and quark- 
antiquark bound states, we can use the partial conservation of axial current (PCAC) 
in constructing the interpolating currents, 9 /J (s(x)7 /i 7 5 M(x)), 

d^{s(x)^^u{x)) = (m s + m u )s(x)i^ 5 u(x), 
(0|^(s(0) 7 M 75«(0))|K(g)) = (m s + m u )(0\s(0)i^ 5 u(0)\K(q)), 

= fxq 2 = fi<m 2 K . 

If we take the s(x) 7^7511 (a;) as the interpolating current, more care has to be taken 
about the possible contaminations from the axial-vector mesons, furthermore, the 
calculation will be more tedious ( with d^{s(x)^^u(x)) ). The matrix element of 
the pseudoscalar current between the vacuum and K state can be taken as 

(0\s(0)ils<0)\K(q)) = X K = f *f K , (5) 

m u + m s 

the values of the depend on the masses of the s and u quarks which have uncer- 
tainties, our numerical results indicate small variations of those masses will not lead 
to large changes about the values of the coupling constant gewx- For the neutron, 
we take the Ioffe current jUj- The Fierz re-ordering of the interpolating currents t]q 
and t]n can lead to the following sub-structures, 

VN = e abc {(u T a Cd b ) l5 d c - (u T a C l5 d b )d c ] , (6) 

Ve = -^^ £abC ( u a C 'y5 d b) {-d c (si^f 5 u) - 7 M 4(^757^) 

-\o^d c (sia^u) - 7^754(^7^) - l<jd c (siu) - (u d) 1 . (7) 
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A naive result of the Fierz re-ordering may be the appearance of the reducible 
contributions with the sub-structure of udd — us (i.e. N — K) clusters in the two- 
point correlation function [TUj . 



however, in our calculations with the interpolating current t]q in Eq.(2), no such 
factorable udd — us terms appear, so there are no reducible N — K contributions 
to the correlation function IT(j>). The re-ordering in the Dirac spin space is always 
accompanied with the color re- arrangement, which involves the underlying dynam- 
ics. If we want to factorize out some N — K contributions from the II (p) ( with 
the current i]® in Eq.(2) ), tedious manipulations about the re-ordering in the color 
and Dirac spin space must be done. There are no direct iV — K ( or udd — us ) 
components in the interpolating current rye, which can readily decay to the NK 
final state and result in large width. If the + (154O) is really a pentaquark state 
not a N — K molecule, as the + (154O) lies above the NK threshold and no need for 
additional quark- ant iquark pairs creation in decay, the decay must be OZI super- 
allowed and the width is supposed be large, say, about several hundred MeV; to 
produce the narrow width, some huge energy barriers are needed to stabilize the 
+ (154O) in case the kinematical interpretation can not work here. The appearance 
of the N — K component in the Fierz re-ordering maybe manifest the possibility ( 
not the probability ) of the evolution from the + (154O) to the NK final state with- 
out net quark-antiquark pairs creation ( maybe the quark-antiquark pairs created 
and annihilated subsequently ), which is significantly in contrast to the conventional 
baryons, however, we have no knowledge about the detailed process of the evolution. 
The G + (1540) may evolve to the N — K final state, and the N — K final state is 
not presented in the components of the initial pentaquark state + (154O), how to 
implement the evolution with small probability involves complex quark-gluon inter- 
actions, whether just re-arrangement in the color space, or creation and annihilation 
of quark-antiquark pairs. In additive constituent quark models, whether or not ad- 
ditional relative P wave is introduced to changed the ground states from negative 
parity to positivity parity, special configurations are needed to take into account the 
narrow decay width as results of small overlaps of the internal and external N — K 
wave- functions [TT], the kinematical interpretations based on the color-flavor-spin ( 
i.e. SU(3) C x SU(3)f x SU(2) S ) group theory resort to the possible small overlaps, 
how to realized the small overlaps needs complex re-ordering in the color-flavor-spin 
space, if there are really no energy barriers to prevent the re- arrangement. While in 
the cluster quark models, typically, the diquark-diquark-antiquark model and 
diquark-triquark model ^3], extra barriers, for example, relative P waves, are intro- 
duced dynamically to prevent the ready decay. In fact, the re-ordering in the color 
and Dirac spin space involves complex strong interactions, and we know little now 
about the dynamics which determine the under-structures of the exotic pentaquark 
states. The mismatches between the color-flavor-spin states in the initial pentaquark 
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and final baryon-meson color singlet can result in suppression of the decay naturally 
[T"H I14j . Due to the spontaneous breaking of the chiral symmetry and the Gold- 
stone nature of the pseudoscalar mesons tt, K and r/, the quarks may have direct 
interactions with the pseudoscalar mesons, which lead to the success of some chiral 
quark models. The dominating interactions which determine the exotic pentaquark 
states be color-spin type A£ • \ C jO % ■ cr 3 or flavor-spin type \{ ■ Ajcr* • are still in 
hot debates [T^], the naive Fierz re-ordering can lead to direct N — K component 
in the B + (1540) will not work here. If there are really some N — K components 
in the interpolating current tjq, they should be factorized out, the remainder can 
not have the correct quantum numbers to interpolate the B + (1540). In the QCD 
sum rules, we construct the interpolating currents with the same quantum number 
as the corresponding mesons and baryons, that is enough; the knowledge about the 
structures of the hadrons can be of much help in the constructing. 

In Ref . jHj , the narrow decay width is attributed to the minor breaking of chirality 
conservation, the color re- arrangement due to the hard gluon-exchange can result 
in strong suppression of the decay, re ~ a 2 < 0|gg|0 > 2 . In the article, we take 
quantitative analysis of the decay width within the framework of the QCD sum rules 
approach. 

The diquark-triquark type interpolating current t)q(x) is more likely related to 
a negative parity pentaquark state, in this work, we make assumption that the 
parity of the + (154O) to be negative and study the decay width with the following 
Lagrangian density, 



C 



igeNK®KN. 



(9) 



According to the basic assumption of current-hadron duality in the QCD sum rules 
approach j3] , we insert a complete series of intermediate states satisfying the unitar- 
ity principle with the same quantum numbers as the interpolating currents jx 
and 7]q into the correlation function in Eq.(l) to obtain the hadronic representation. 
After isolating the double-pole and single-pole terms of the lowest ground states, we 
get the following result [SJ US] , 



T(p,q) 



\ ttIq — p' 2 {m 2 N — p 2 ) {m 2 K — q 2 ) 
W^q^Pu +•••} + •••, 
f 9&nk A© A jv 



+ 



A{p' 2 ,g 2 ) B(p' 2 ,p 2 ) 

1 o 1 < 



ni 



N 



m 



K 



rn 



e 



P 



/2 











p 2 )(t — q 2 ) 



s 

00 dc A(p' 2 ,q 2 )5(s-m 2 N ) | ^ dt B(p' 2 lP 2 )5(t-m 2 K ) 



s-p 
Wq»Pv +•••} + • 



t- q 2 



+ 



+ 



(10) 
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with 



WV) 



00 dt^^ 
,2 



00 



ds 



t-q 2 ' 

Pb(p' 2 ,s) 
s — p z 



(12) 



where the following definitions have been used, 



(0\ VN \N(p,s)) 
(e(p',s')\fj e \0) 



X N u(p, s) 
X @ u(p', s') 



The coupling constants Xn and Ae can be determined from the two-point QCD sum 
rules, for the Ae, we use the correlation function H(p), substitute the tjn for the rje in 
Eq.(8), we can obtain the A^. In this article, we choose the Dirac tensor structure 
<7 MI ' q^Pv for analysis, while in Ref.f.6 a , the authors take the structure js&^q^Pv ■ 
Here the residues of the single-pole terms A(p' 2 ,q 2 ) and B(p' 2 ,p 2 ) have complex 
dependence on the transitions between the ground states and high resonances ( or 
continuum states ). We have no knowledge about the transitions, even the existence 
of the B + (1540) is not firmly established. However, the contributions from the pole- 
continuum transitions are not exponentially suppressed compared with the double- 
pole terms, even after double Borel transform, furthermore, the contributions can 
be as large as or larger than the double-pole terms and must be explicitly included 
in the sum rules. We only have the fact that the B + (1540) lies a little above the 
NK threshold, the contributions from the + (154O) can be factorized out, so the 
spectral densities pa and pb can be parameterized as 



pA(p' 2 ,t) 

Pb{p' 2 ,s) 



EE(p' 2 ,t) 

2 /9 ' 

m Q+ ~~ P 

FF(p' 2 ,s) 

m Q+ ~ P' 2 



(13) 
(14) 



The two unknown functions EE and FF have complex dependence on the transitions 
between the ground states and high resonances ( or continuum states ). From the 
Eqs.( 10-14), we can obtain 



T{p, q) 



) 90NkXbXnXk 



1 



V 



m e ~ P' 2 



+ 



m % ~ P' 2 



{m 2 N -p 2 ){m 2 K - q 2 ) 
EE(p' 2 ,t) 



m 2 N — p 2 



dt 



m 2 K — q 2 



p A 



Wq.Pv +•••} + 
) geNKXeXxXK 



\ m e — p' 2 {m 2 N — p 2 ) (m 
W" u q,Pu +■■■} + 



k - Q 2 : 



+ 



1 



cc 



DD 



m e ~ V Y m N ~ V m K ~ 1 



here we introduce two constants CC and DD for convenience, 



CC 



= r as) 



Taking the CC and DD as some unknown constants has smeared the complex de- 
pendence on the energy and high resonance masses ( or continuum states ), which 
will certainly impair the prediction power. We have no knowledge about the tran- 
sitions between the pentaquark state G + (1540) and the excited states ( or high 
resonances, or continuum states ), the CC and DD can be taken as free parameters 
, we choose the suitable values for the CC and DD to eliminate the contaminations 
from the single-pole terms to obtain the reliable sum rules. The contributions from 
the single-pole terms may as large as or larger than the double-pole term, in practical 
manipulations, the CC and DD can be fitted to give stable sum rules with respect 
to variations of the Borel parameter M 2 in a suitable interval. If the final numerical 
results are insensitive to the threshold parameters s , t and there really exists a 
platform with the variations of the Borel parameters M 2 and Mf , the predictions 
make sense. 

The calculation of operator product expansion in the deep Euclidean space-time 
region is straightforward and tedious, here technical details are neglected for sim- 
plicity, once the analytical results are obtained, then we can express the correlation 
functions at the level of quark-gluon degrees of freedom into the following form 
through dispersion relation, 

We choose p 2 = p' 2 = —P 2 and q 2 = —Q 2 , then take double Borel transform with 
respect to the variables P 2 and Q 2 respectively, match Eq.(16) with Eq.(19), finally 
we obtain the sum rules for the strong coupling constant gewK, 

2 _ m l _ m N 2 2 

_ m K g Mf _ g Mf _ m K _ m C-) 

^e^A^AjvAee ^ = 5— + Ce ^ e ^ = 



rria — m 



N 



^f 21m s MfM 2 E 2 (s)E (t) llm s (qq) 2 M 2 E (t) 
\ 2 n 4!7r 6 2 7 3tt 4 

7 l(qq) + (55)] M?E 2 (s) 11 [(qq) 3 + (qq) 2 (ss)] 



2 8 4!tt 6 2 4 3% 4 



(20) 
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where 




Here the C ( proportional to the DD, as the CC terms are eliminated ) denotes 
the contributions from the transitions between the ground and excited states ( or 
high resonances, or continuum states ), we can choose the suitable values for C to 
eliminate the contaminations to obtain the stable sum rules with the variations of 
the Borel parameters Mf and Mf . 



3 Numerical Results 

The parameters for the condensates are chosen to be the standard values 4J, (ss) = 
(0.8 ± 0.1)(uu), (qq) = (uu) = {dd) = -(240 ± lOMeV) 3 , m u = m d = and 
m s = (140 ± 10)MeV. Small variations of those condensates will not lead to large 
changes about the numerical values. The coupling constants are taken as = 
(2.4 ± 0.2) x 10~ 2 GeV 3 HQ EE3 and A© = (1.4 ± 0.2) x 10~ 4 GeV 6 [7 from the 
two-point QCD sum rules, see Eq.(8) for example. The threshold parameters sq 
and to are chosen to vary between (1.8 — 2.0)GeV 2 and (0.8 — 1.0)GeV 2 respectively 
to avoid possible contaminations from high resonances and continuum states. The 
Borel parameters are taken as Mf = (1.0 - 1.5)GeV 2 and Mf = (1.0 - 2.0)GeV 2 to 
obtain the stable sum rules. Finally we obtain the values for the coupling constant 

\g<dNK\i 

\g eNK \ =0.175 ±0.084, (21) 



r@ = S^T^entfKraJv + m e ) 2 -m 2 K ]^ A(m|, m%, m 2 K ) , 



< AMeV, (22) 

w2 2 2\ /2, 2 2\2 a 2 2 

A(m e ,m N ,m K ) = {m e + m N -m K ) - Am e m N . 

which can explain the narrow width V < lOMeV naturally. The values of the 
coupling constant geNK with the variations of the threshold parameters (so, t ) 
and Borel parameters (Mf , Mf) are shown in Fig.l and Fig. 2 for (ss) = 0.8{uu), 
(qq) = -(240Me1/) 3 , m s = UOMeV. 



4 Conclusion 

In this article, we take the point of view that the pentaquark state G + (1540) has 
negative parity, and choose the diquark-triquark type interpolating current to cal- 
culate the strong coupling constant geNK within the framework of the QCD sum 
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s =1.8GeV 2 ; 
s =1.9GeV 2 ; 
s =2.0GeV 2 . 



1.2 1.4 1.6 1.8 2.0 

M^GeV 2 ) 

Figure 1: \g eNK \ with for M| = 1.3GW 2 , t = 0.9Ge^ 2 . 



t =0.8GeV 2 ; 
t =0.9GeV 2 ; 
t =1.0GeV 2 . 



1.1 1.2 1.3 1.4 1.5 

M 2 2 (GeV 2 ) 

Figure 2: | 5eAri r| with M| for M\ = UGeV 2 , s = 1.9GeV 2 
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rules approach. Our numerical results indicate that values of the strong coupling 
constant geNK are ver Y small, \geNx\ = 0.175 ± 0.084, and the width r < 4MeV, 
which can explain the narrow width T < lOMeV naturally. 
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